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The static and dynamic properties of V*^ spins {S — 1/2) in the frustrated square lattice com- 
pound Pb2(VO)(P04)2 were investigated by means of magnetic susceptibility x and "^^P nuclear 
magnetic resonance (NMR) shift (K) and ''^P nuclear spin-lattice relaxation rate l/Ti measurements 
on a single crystal. This compound exhibits long-range antiferromagnetic order below Tn — 3.65 K. 
NMR spectra above Tjv show two distinct lines corresponding to two inequivalent P-sites present in 
the crystal structure. The observed asymmetry in hyperfine coupling constant for the in- plane (PI) 
P-site directly points towards a distortion in the square lattice at the microscopic level, consistent 
with the monoclinic crystal structure. The nearest- and next-nearest-neighbor exchange couplings 
were estimated by fitting K versus temperature T by a high temperature series expansion for the 
spin susceptibility of the frustrated square lattice to be Ji/fes = ( — 5.4±0.5) K (ferromagnetic) and 
J2/kB ~ (9.3 ± 0.6) K (antiferromagnetic), respectively. l/{TiTx) is almost T-independent at high 
temperatures due to random fluctuation of spin moments. Below 20 K, the compound shows an 
enhancement of l/{TiTx) which arises from a growth of antiferromagnetic spin correlations above 
Tn. Below Tn and for the field applied along the c-axis, the NMR spectrum for the PI site splits into 
two satellites and the spacing between them increases monotonically with decreasing T which is a 
direct evidence of a columnar antiferromagnetic ordering with spins lying in the a6-plane. This type 
of magnetic ordering is consistent with expectation from the J2/J1 — —1.72 ratio. The critical expo- 
nent (3 = 0.25 ±0.02 estimated from the temperature dependence of the sublattice magnetization as 
measured by ^^P NMR at 11.13 MHz is close to the value (0.231) predicted for the two-dimensional 
XY model. 

PACS numbers: 75.50.Ee, 75.40.Cx, 71.20.Ps 



I. INTRODUCTION 

Understanding the ground state properties of frus- 
trated low-dimensional spin systems is a central issue in 
current condensed matter physics. These systems have 
competing interactions which stabilize different states 
with distinct symmetries.^ An interesting phenomenon 
in spin systems is the formation of a spin-liquid with no 
long-range magnetic order (LRO) due to suppression by 
geometric/magnetic frustration. A simple example is the 
frustrated S = 1/2 square lattice (FSL) model. In this 
model which is also known as Ji — J2 model, the spin 
Hamiltonian is 

n = (fee Ji) ^ ■ Sj + (M) J2 • ^"fc' 

<ij> <ik> 

where the first sum is over nearest-neighbor spin pairs 
and the second is over next-nearest- neighbor spin pairs. 
A positive J corresponds to an antiferromagnetic (AF) 
exchange interaction and a negative J to a ferromagnetic 
one. In this paper, Ji and J2 are expressed in tempera- 
ture units, and fee is Boltzmann's constant. The nearest- 
neighbor (NN) interaction Ji along the side of the square 
can compete with the next-nearest-neighbor (NNN) in- 
teraction J2 along the diagonal of the square if, e.g., both 
are AF. Recently a generalized phase diagram including 
both antiferromagnetic and ferromagnetic Ji and J2 has 



been proposed theoretically based on the frustration ra- 
tio a = J2/J1 or frustration angle = tan~^(J2/Ji) 
as shown in Fig. [TJMsi It has three different ordered 
phases: ferromagnet [FM, wave vector Q — (0,0)], Neel 
antiferromagnet [NAF, Q = (7r,7r)], and columnar an- 
tiferromagnet [CAE, Q = (7r,0) or (0,7r)]. The critical 
regions with disordered ground states are predicted to 
occur at the boundaries NAF-CAF and the CAE-EM 
which are quantum spin liquid (QSL) phases. A gap- 
less nematic state is suggested for a ^ —0.5,"^ while dif- 
ferent dimer phases (including resonating-valence-bond- 
type ones) are claimed to exist for a close to CSj^iii^ 
In addition to its intrinsic interest, a better understand- 
ing of the Ji — J2 model is expected to play a vital role 
in understanding the magnetism in recently discovered 
FeAs-superconducting parent comp ounds 1^ 1 ^ ° 1 ^ ^ 

Experimentally, however, few compounds satisfying 
the Ji — J2 square spin lattice model have been realized 
so far. Li2VOA:04 {X = Si, Ge}^^ are two well-studied 
compounds with antiferromagnetic Ji and J2 and lie in 
the CAE state of the phase diagram, far away from the 
QSL regime. ^"^'^^ Recently two new perovskite type com- 
pounds (CuAr)LaNb207 {X — CI, Br) were claimed to 
realize the Ji — J2 model with ferromagnetic Ji and an- 
tiferromagnetic J2'^'^ However the validity of the frus- 
trated square lattice model for these two systems is still 
unclear Recent band structure calculations suggested 
that despite the layered structure of (CuCl)LaNb207, 
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FIG. 1: (Color online) The Ji — J2 phase diagram of the 
two-dimensional spin- 1/2 square spin lattice showing different 
ordered phases and the quantum spin liquid (QSL) regimes. 
The two QSL regimes are the two vertically hatched areas. 
The boundaries of the QSL regimes are not precisely known. 
The other regions of the phase diagram contain the ferro- 
magnetic (FM), columnar antiferromagnetic (CAF), and Neel 
antiferromagnetic (NAF) phases. The location of all the in- 
vestigated compounds so far in the columnar antiferromag- 
netic (CAF) phase are shown.— i^^i^i^ The radius of the cir- 
cle is 12.2 K. Due to relatively large exchange couplings, 
PbVOa and VOM0O4 are not included in the phase diagram, 
for which Ji = 203 K and J2 = 77 K (Ref.|2l}and Ji = 110 K 
and J2 = 22 respectively. (After Refslii andlH) 



the spin system has rather pronounced one-dimensional 
character .ia The compounds VOM0O4 (Refs. HO andllH) 
and PbVOa (Ref. HI) have been studied in light of the 
Ji — J2 model and were deduced to lie in the NAF region 
of the phase diagram. 

The vanadium phosphates AA' VO(P04)2 {AA' = Pb2, 
SrZn, BaZn, and BaCd) present another realization of the 
Ji — J2 model with ferromagnetic Ji and antiferromag- 
netic J2 and fall in the CAF region of the phase diagram 
close to the QSL regime i^i^iiS^i^^ The magnetic proper- 
ties of these compounds have been recently investigated 
by means of magnetization,^>2& heat capacity ) ^^1^^ neu- 
tron scattering , ^^'^^ and muon spin resonance^S (/iSR) 
experiments. All these measurements were done on poly- 
crystalline materials. However a clear picture about 
the spin dynamics in these compounds requires measure- 
ments on single crystals. Here we present '^^P NMR and 
magnetic susceptibility x measurements on a large single 
crystal of Pb2VO(P04)2- 

Pb2(VO)(P04)2 crystallizes in a monoclinic struc- 
ture with space group P2i/a and lattice parameters 
a = 8.747(4) A, & = 9.016(5) A, c = 9.863(9) A, and 
P = 100.96(4)°,^ The structure of Pb2(VO)(P04)2, pre- 
sented in the upper panel of Fig. [21 is composed of 
[VOPO4] layers extending parallel to the a6-plane which 
are separated along the c-axis by Pb atoms and isolated 
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FIG. 2: (Color online) Upper panel: schematic diagram of 
[VOPO4] layers in Pb2(VO)(P04)2 running parallel to the 
afo-plane. Lower panel: projection onto the ab-plane of the 
unit cell showing the arrangement of VO5 pyramids and PIO4 
tetrahedra in the plane. The approximate square lattice 
formed by four V^^ spins (labeled as A, B, C, and D) is also 
shown indicating the nearest-neighbor (Ji) and next-nearest- 
neighbor (J2) exchange interactions. 



PO4 tetrahedra as shown. The layers are slightly modu- 
lated due to lack of tetragonal symmetry. It is important 
to mention here that there are two inequivalent phos- 
phorus sites (PI and P2) present in the crystal structure. 
The [VOPO4] layers are formed by corner-sharing V^+Os 
pyramids with PIO4 tetrahedra while the other site forms 
isolated P2O4 tetrahedra lying between the [VOPO4] lay- 
ers. Thus the PI site is expected to be strongly coupled 
and the P2 site weakly coupled to the (3d\ 5 = 1/2) 
spins. In the lower panel of Fig. [2] a section of the ab- 
planc is shown. In the approximate V4 square A-B-C-D 
shown, the V*"^ spins interact via a V-O-Pl-O-V path- 
way involving two inequivalent oxygen atoms. The dis- 
tances and angles of the pathways between the different 
V spins are also different which reflects the inequivalent 
superexchange paths. We have listed in Table U the dis- 
tances between the pairs of four spins A, B, C, D, 
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TABLE I: The distance between pairs of V spins along the 
edges and along the diagonals, and the interior angle between 
them, for the approximately square lattice shown in the lower 
panel of Fig.gjfor Pb2VO(P04)2. 



V-V pair 


AB 


BC 


CD 


DA 


AC BD 


Edges (A) 


4.669 


4.665 


4.669 


4.410 




Diagonals (A) 










6.267 6.295 




/ABC 


/BCD 


/CD A 


/DAB 




Angles (°) 


84.34 


84.82 


87.25 


87.75 





and the angles between them. Small differences are seen 
in the distances along the edges, the angles between the 
edges, and the distances along the diagonals from square 
symmetry, indicating a small distortion of the square. 

The exchange interactions Ji and J2 between V spins- 
1/2 in Pb2VO(P04)2 were reported to be -6.0 K and 
9.8 K, respectively, and the compound was reported to 
undergo antiferromagnetic ordering at Tn — 3.65 K, 
likely of CAF typei^ Recently Tsirlin and Rosner— 
pointed out theoretically from band structure calcula- 
tions that due to the reduced symmetry the Ji and J2 
values within the V4 square are multivalued: (Jj^, Jj^ ) 
and ( J2, J2 )^ respectively and the experimental estimates 
of Ji and J2 should be considered as averaged NN and 
NNN couplings: Ji = {j[ + jI)/2 and J2 = [j'2 + J2 )/2- 
From neutron scattering experiments, a diffuse magnetic 
background was observed in the paramagnetic state at 
T = 20 K which persists well below Tn- The ground 
state was suggested to be CAF type with a reduced or- 
dered moment (/x ~ 0.5 /is/V"'"), where /xb is the Bohr 
magneton . ^^'^^ For a spin- 1/2 with (7-factor g — 2, one 
would expect an ordered moment fi = gSfiB = 1 Mb per 
V atom. The diffuse background and the reduced or- 
dered moment may be indications of quantum disorder 
predicted theoretically near the QSL regime.'^ /iSR ex- 
periments also show a broad distribution of field at the 
muon site below Tn suggesting possible structural and/or 
magnetic disorder. The temperature dependence of the 
longitudinal muon relaxation rate above Tn and of the 
order parameter below Tn agree with expectations for 
the 2D XY model. Both neutron scattering and /iSR 
experiments were performed on powder samples. In un- 
published work, the x(T) of a single crystal above Tn 
was found to be isotropic along different orientations. 
Below Tn , it is strongly anisotropic with the 6-axis being 
the easy axis. ^4 

Nuclear magnetic resonance (NMR) is a very powerful 
experimental tool to extract microscopic information at 
the individual sites of a crystal. The technique has been 
applied to many transition metal oxides and has played a 
significant role in elucidating their microscopic magnetic 
characters. In this paper, we report a detailed exper- 
imental investigation of static and dynamic properties 
of spins in Pb2(VO)(P04)2 by means of ^ip NMR 
and X measurements. We have carried out '^^P-NMR not 



only on a polycrystalline sample but also on a large single 
crystal in order to shed light on the spin structure and 
possible quantum disorder in Pb2(VO)(P04)2- 

II. EXPERIMENTAL 

Syntheses of the polycrystalline sample and of the crys- 
tal are reported in Refs. [2^ and [s^, respectively. The 
NMR measurements were carried out using pulsed NMR 
techniques on "^^P (nuclear spin I — 1/2 and gyro- 
magnetic ratio 7jv/27r — 17.237 MHz/T) nuclei in the 
temperature range 1.5 K < T < 300 K. Lane x-ray 
back-scattering measurements previously determined in 
Ref. [U the orientation of the crystal axes with respect to 
the crystal faces for the crystal on which NMR measure- 
ments were performed. We have done the NMR measure- 
ments at two different radio frequencies of 70 MHz and 
11.13 MHz which correspond to an applied field of about 
4.06 T and 0.65 T, respectively. Spectra were obtained 
either by Fourier transform of the NMR echo signals or by 
sweeping the field. The NMR shift K = {1/ — v^cf)/vref 
was determined by measuring the resonance frequency 
of the sample {v) with respect to nonmagnetic reference 
H3PO4 (resonance frequency J^rof)- The ^^P spin-lattice 
relaxation rate 1/Ti was measured by the conventional 
single saturation pulse method. For the analysis of NMR 
data, we measured the magnetic susceptibility x{T) of 
the single crystal at 4 T for field applied along c-axis in 
a commercial (Quantum Design) SQUID (Superconduct- 
ing Quantum Interference Device) magnetometer. 

III. RESULTS 
A. NMR spectra above Tn 

The '^^P nucleus has nuclear spin / — 1/2. For a / = 
1/2 nucleus one expects a single spectral line for each 
inequivalent P sitei^ As shown in the crystal structure 
(Fig. [21), Pb2VO(P04)2 has two inequivalent P sites. We 
indeed observed two narrow spectral lines corresponding 
to the two P sites for all three orientations of the crystal 
above Tn (e.g., see ^ip NMR spectra at 4.27 K in Fig.©. 
For each P site the line position was found to shift with 
temperature. However for the strongly coupled PI site 
the shift was found to be much stronger compared to the 
weakly coupled P2 site and is also strongly orientation 
dependent. In the present work we have mainly focused 
on the PI site and therefore data for the P2 site are not 
shown. Figure [3] shows the "^^P NMR spectra for the PI 
site measured on a single crystal at different temperatures 
and field applied along different orientations. 

FigureSKa) presents the temperature dependence oi K 
for the PI site derived from the data in Fig. [3l It shows 
a strong anisotropy along different directions. Along the 
a- and &-axes it shows a strong positive shift while along 
the c-axis it shows a strong negative shift. On the other 



4 




70.0 70.2 70.4 70.6 70.8 



V (MHz) 



FIG. 3: (Color online) High-field Fourier transform ^^P 
NMR spectra for the PI site at different temperatures T 
(T > T-m) and along different orientations for single crystalline 
Pb2(VO)(P04)2. The vertical dotted line corresponds to the 
^^P resonance frequency of the reference sample H3PO4. 
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FIG. 4: (Color online) (a) Temperature dependent NMR shift 
K vs. T for the PI site along different field orientations. The 
solid line is the fit of Ki^o by Eq. (b) ^^P shift K vs. x 
measured at 4 T is plotted with temperature as an implicit 
parameter for all three orientations. The solid lines are linear 
fits. 



hand for the P2 site the shift is very weak. The stronger 
and weaker shifts for the PI and P2 sites, respectively, 
suggest that the former one is strongly coupled while the 
latter one is weakly coupled to the V''+ spins, consistent 
with expectation from the crystal structure. The temper- 
ature dependence of the isotropic ^^^Pl NMR shift Ki^o 
was calculated as i^iso = {Ka + + Kc)/3 and is also 
shown in Fig. [D^a). ft shows a marked T dependence 
like x{T)-'^'^ At high temperatures Ki^o varies in a Curie- 
Weiss manner and then passes through a broad maximum 
at around 9 K reflecting low-dimensional short-range an- 
tiferromagnetic ordering. Since x(^) for the single crys- 
tal at T > Tn is isotropic along different orientations;^ 
this anisotropic K(T) is likely due to asymmetry in the 
hyperfine coupling constant between the V*"*" spins and 
the P nuclear spins. 

A ^^P spectrum on a polycrystalline sample is shown in 
Fig. [5l Along with the most intense central line there is a 
broad background containing extra shoulder-like features 
on either side. In an attempt to fit the experimental spec- 
tra we assumed that the central part with a small asym- 
metry comes from the weakly coupled P2 site and the 
broader one with large asymmetry is from the strongly 



coupled PI site. In this way the superposition of these 
two spectra gives a reliable fit as shown in Fig. [51 One 
can see in Fig. [5] that for the broad spectrum there are 
three shoulder-like features clearly visible. These shoul- 
der positions were found to shift with temperature. We 
took the derivative of the spectra and picked the point 
of inflection for each shoulder as a function of tempera- 
ture (not shown) and they overlap nicely with the single 
crystal K{T) confirming that the shoulders are due to 
anisotropic shift at the PI site. On the other hand the 
narrow line (P2 site) was found to shift very weakly with 
temperature. 

The NMR shift is a direct measure of the spin suscepti- 
bility Xspin, and quite generally K(T) is written in terms 

of xspi„(r) nsM 

K{T)^Ko + ^Xspin{T), (2) 

where Kq is the temperature-independent chemical (or- 
bital) shift, Na is the Avogadro number, ^hf is the hy- 
perfine coupling constant between the P nuclear spins 
and the V^"*" electronic spins, and Xspin/A^A is expressed 
in units of (Ib/Oc per electronic spin, which in our case 
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FIG. 5: (Color online) The ^^P NMR spectrum at 14.8 K mea- 
sured at 70.1218 MHz on a polycrystalline Pb2(VO)(P04)2 
sample (open circles). The dashed and dotted lines are the 
theoretical fits for the PI and P2 sites, respectively, and 
the solid line is the total fit which is a superposition of 
the PI and P2 fits. The NMR shift values along the a- 
, h-, and c-directions obtained from the fitting are [Ka — 
5500 ppm, Kb ~ 3430 ppm, and Kc ~ —4160 ppm) and 
{Ka — 1960 ppm, Kb ~ —342 ppm, and Kc — —512 ppm) for 
the PI and P2 sites, respectively. 



means per formula unit. The conventional scheme for 
calculating Am is to take the slope of the K vs Xspin 
plot with T as an implicit parameter. Here we mainly 
focus on the PI site. For each orientation the K vs. 
X plot is fitted very well by a straight line [Fig. [IJb)] 
over the whole temperature range (T > Tn) yielding 
Am = (1882 ±40), (1251 ±42), and -(1642 ±55) Oe/^B 
along the a-, 5-, and c-directions, respectively. The 
isotropic hyperfine coupling constant is then calculated 
to be Aso = {Aa + Ab + ^e)/3 = (497 ± 46) Oc/mb- 

In order to estimate the exchange couplings Ji and J2, 
we fitted the temperature dependence of Kiso above 15 K 
by Eq. ^ , where Xspin is the the high-temperature series 
expansion prediction for the molar spin susceptibility of 
the frustrated square lattice model given by^^ 



Xspin (^) 



knT ^ 



knT 



where Cm,n are the coefficients listed in Table I of Ref. llSl. 
g is the (7-factor, and is the Boltzmann constant. 
From the electron paramagnetic resonance measurements 
at room temperature, g parallel {g\\) and perpendicu- 
lar {g±) to the c-axis were found to be 1.929 and 1.966, 

respectively!^ The isotropic g = 



'II I 2g^)/3 was cal- 
culated to be 1.95 which was kept fixed during the fitting 
procedure. In this way we obtained Kq = (462±11) ppm. 
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FIG. 6: (Color online) (a) Nuclear spin-lattice relaxation rate 
1/Ti vs. temperature T measured for the PI site at 70 and 
11.13 MHz with the magnetic field applied along the c-axis. 
Below Tn measurements were done at the central field corre- 
sponding to the right-hand-side satellite [see Fig. [Tjb)] . The 
solid and dashed lines represent T''' and behaviors, respec- 
tively, (b) l/(xrir) (left y-axis) and x (right y-axis) plotted 
as a function of T . The inset shows the H vs. T phase dia- 
gram where Tn values obtained from our l/Ti measurements 
at the PI site and the values from heat capacity Cp measure- 
ments (Ref. lU) are plotted. 



Am = (481 ± 20) Oe/^iB, -h = (-5.4 ± 0.5) K, and 
J2 = (9.3 ±0.6) K. The fit is shown in Fig.HJa) as a solid 
line. These values of J\ and J 2 are close to the previously 
reported values estimated from x(T) analysis^^iSl and 
also consistent with the saturation field obtained from 
high-field magnetization isotherm measurements.^^ The 
estimated ^hf matches with the ^iso calculated above 
from the K-x analysis within the error bars. 



B. Nuclear spin-lattice relaxation rate l/Ti 

The "^^P nuclear spin-lattice relaxation rate 1/Ti was 
measured with the magnetic field applied parallel to the 
c-axis. For a / = 1/2 nucleus the recovery of the lon- 
gitudinal magnetization is expected to follow a single- 
exponential behavior. In Pb2VO(P04)2, the recovery of 
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the nuclear magnetization after a saturation pulse was 
indeed fitted well by the exponential function 
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1 - 



Mn 



(4) 



where M{t) is the nuclear magnetization at a time t af- 
ter the saturation pulse and Mq is the equilibrium nu- 
clear magnetization. The temperature dependence of 
1/Ti measured for the PI site for H \\ c is presented 
in Fig. [eja). At high temperatures (T > 20 K), 1/Ti 
is temperature independent. In the high temperature 
limit T ^ J, a temperature- independent 1/Ti behav- 
ior is typical when paramagnetic moments fluctuate fast 
and at random.— With decrease in temperature, 1/Ti 
decreases slowly for T < 15 K and then shows a peak 
around 3.93 K. This decrease is very similar to that 
observed previously in the cases of the antiferromag- 
netic square lattices VOM0O4 from ^^Mo NMR^o and 
[Cu(HC02)2-4D20] from NMR where the decrease of 
1/Ti above Tn is explained by cancellation of the an- 
tiferromagnetic spin fluctuations at the positions of the 
probed nuclei.'^- Below the peak 1/Ti again decreases 
smoothly towards zero. 

For the polycrystalline sample, l/Ti measurements 
were not possible due to broad spectra and overlapping 
of signals from the two P sites. 



C. NMR spectra below Tn 

Below Tn both the PI and P2 lines were found to split 
(Fig. [7|) indicating that both P sites are experiencing the 
static internal field in the ordered state. The PI site 
line splits into two satellite lines while the P2 site line 
splits into four satellite lines as a result of the hyperfine 
field between the P nuclei and the ordered moments. 
For the 70 MHz measurements [Fig. [Tlja)], the left-hand- 
side satellite of the PI site overlaps with the four split 
satellites of the P2 site making it difficult to distinguish 
the peaks. It is possible to separately measure the spectra 
for the two different sites by measuring at low frequency 
and using a faster repetition rate since the relaxation 
rates associated with the two sites are different (for the 
PI site Ti is shorter than for the P2 site). Therefore we 
remeasured the spectra for the PI site at 11.13 MHz using 
a faster repetition rate and the results are presented in 
Fig.Hb). 

The internal field -ffint, which is proportional to the 
V^^ sublattice magnetization, was determined from half 
the separation between the split lines. The temperature 
dependences of i7int are shown in Fig. [5] for both the P- 
sites measured for H \\ c. For the PI site at 70 MHz, the 
left-hand-side satellite was chosen as the peak of equal 
intensity to the right-hand-side satellite. The Hi^tiT) for 
the PI site was found to be almost independent of field 
below Tn. In order to extract the critical exponent of the 
order parameter (sublattice magnetization), H\i^t{T) was 
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FIG. 7: (Color online) Temperature dependent ^^P NMR 
spectra measured at (a) 70.00 MHz and (b) 11.13 MHz for 
i7 |j c. The spectra in the paramagnetic state split below Tn 
into two and four lines for the PI and P2 sites, respectively, 
due to the internal field iifint. The dashed lines are guides to 
the eye for the shifts of the satellite lines. 
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FIG. 8: (Color online) Temperature T dependence of the in- 
ternal field i/int obtained from NMR spectra for i7 || c in the 
ordered state (a) for the PI site measured at 11.13 and 70 
MHz and (b) for the P2 site measured at 70 MHz. i7int is 
proportional to the V*^ sublattice magnetization. The solid 
lines are fits by Eq. to the data, as described in the text. 



fitted by the power law 

H,nt{T) ^ Ho (^1 - . (5) 

One can notice tliat Hi^^ decreases sliarply on approacli- 
ing Tn- For an accurate determination of the critical 
exponent /?, one needs data points close to T/v (i.e. in 
the critical region). We have estimated /3 by fitting the 
data points in temperature steps of 0.01 K as close as 
possible to Tn as shown in Fig. [S] For the PI site at 
11.13 MHz, the maximum value of f3 = 0.25 ± 0.02 with 
Tn — 3.655 K was obtained by fitting the data points 
(3.58 K to 3.63 K ) close to Tn. By increasing the number 
of fitting points towards low-Ts, the (3 value was found 
to decrease and the minimum value of /3 = 0.15 ± 0.02 
with Tn — 3.64 K was obtained by fitting the whole tem- 
perature range (1.58 K to 3.63 K). Similarly for the PI 
site at 70 MHz, the maximum value oi P — 0.23 ± 0.03 
with Tn — 3.933 K was obtained by fitting data points 
(3.9 K to 3.93 K) close to Tn and it decreases when low- 
T data points were included. At low-Ts, Hmt develops 
the tendency of saturation and it saturates much faster 
than expected from mean field theory [see the deviation 
of fits in Fig. [5]Ja) at low-Ts]. On the other hand for the 
P2 site at 70 MHz, close to Tn, due to overlapping of the 
lines from the two P sites an accurate determination of 
Hint was not possible; therefore the data points shown in 
Fig. El^b) for the P2 site have large error bars. 

For the polycrystalline sample below T/y, the NMR line 
at 70 MHz was found to broaden abruptly. In order to 
check whether any extra features could be resolved, we 
remeasured the spectra at a low frequency (12 MHz) and 
again with a fast repetition rate. It was found that above 
T/v the line remains narrow and immediately below T/y 
it starts broadening (see Fig. [9]) . The broad line arising 
from the PI site develops two satellites on either side of 
the central peak that are due to the P2 site. With a 
decrease in temperature these two satellites move away 
from each other.— We then plotted half the distance be- 
tween the two satellite positions (-ffint) as a function of 
temperature (inset of Fig. [5]). Fitting the data from 1.7 K 
to 3.6 K by Eq. © yielded (3 = 0.25 ± 0.07. 



IV. DISCUSSION 
A. Statics 

The exchange couplings extracted from our NMR 
iso-shift data are consistent with the values reported 
from x(T) analysis i^^i^^ Based on the J2/J1 ratio, 
Pb2VO(P04)2 with a ~ -1.72 is located in the CAF 
regime of the phase diagram in Fig. [1] In the crys- 
tal structure, squares are formed via a V-O-Pl-O-V su- 
perexchange interaction path. Since PI is strongly cou- 
pled to the V^^ spins via hyperfine interaction, any 
change in the structure properties would be reflected in 
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H(T) 

FIG. 9: (Color online) Temperature-dependent ^^P NMR 
spectra measured at 12 MHz on a polycrystalline sample be- 
low Tn. The dotted lines indicate the splitting of the PI line. 
Inset: Internal field Hint as a function of temperature T. The 
solid line is a fit by Eq. ((5]). 



the '^^P NMR. The total hyperfine coupling constant at 
the PI site is generally the sum of transferred hyperfine 
(^trans) and dipolar (Adip) couplings produced by 
spins, i.e. ^hf = z Atrans + ^dip whcrc z is the num- 
ber of nearest-neighbor V^^ spins of the PI site. The 
anisotropic dipolar couplings were calculated for three 
orientations of the crystal using lattice sums with an as- 
sumption of 1 /zb magntic moment on each V*^ ion, 
which are A'^-^ = 150 Oe//iB, ^dip = 170 Oe//iB, and 
Ajjp — —320 Oe//iB along the a-, 6-, and c-directions, 
respectively. Thus by subtracting the dipolar coupling 
from the total hyperfine coupling, the transferred hyper- 
fine coupling was obtained to be ^trans — 1730, 1080, 
and —1320 Oe//iB, respectively, suggesting that the dom- 
inant contribution to the total hyperfine coupling is due 
to the transferred hyperfine coupling at the PI site. 
The transferred hyperfine coupling at the PI site arises 
mainly from interactions with the four nearest neighbor 
V spins in the plane. The isotropic and anisotropic trans- 
ferred hyperfine couplings originate from P(3s)-0(2p)- 
V(3(i) and P(3p)-0(2p)-V(3d) covalent bonds, respec- 
tively. Since PI is surrounded by four V ions forming an 
approximate square lattice in the plane, the experimen- 
tally observed asymmetry in total hyperfine field in the 
plane indicates inequivalent P(3p)-0(2p)-V(3d) bonds for 
the four nearest-neighbor V ions and hence a distortion 
in the square lattice, consistent with the low symmetry of 
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the crystal structure as pointed out in Sec. I. Our obser- 
vations seem to be consistent with the recent theoretical 
calculation by Tsirlin et alj^ where they have estimated 
the AFM NNN exchange couplings and . The ra- 
tio is J2 / J2 = 0.67, smaller than the value of unity ex- 
pected for a regular frustrated square lattice, consistent 
with a distortion in the square. Theoretical studies have 
indicated that the frustrating next-nearest-neighbor cou- 
pling J2 can lead to a lattice distortion due to spin-lattice 
coupling. ■^^'^^ 

Now we discuss the spin structure in the magnetically 
ordered state. As shown in Fig. [71 for H \\ c, the PI 
line splits into two lines rather than being simply shifted 
to lower field in the antiferromagnetically ordered state. 
We also carried out spectral measurements below Tn for 
the field applied along a- and 6-axes at 11.6 MHz. For 
both directions, the NMR line was also found to split but 
the splitting is much weaker compared to that with the 
field along the c-axis. For instance, at T « 2.9 K, the in- 
ternal field Hint along the a- and 6-axes are ~ 5 Oe and 
^ 10 Oe, respectively, while along the c-axis it is about 
250 Oe.'*'^ Thus the direction of the internal field Hi^t at 
the PI site is almost parallel (or antiparallel) to the ap- 
plied field H along the c-axis. In this case the effective 
field at the PI site would be HcS = H zL Hjat- The satel- 
lite lines are rather sharp with no striking broadening. 
Thus the splitting of the NMR line with no broadening 
is direct evidence of a commensurate magnetic structure. 
If the magnetic structure were incommensurate with the 
lattice, the internal field would be distributed and the 
spectrum would not exhibit sharp resonance lines as seen 
in Fig. [71 As pointed out before by neutron scattering ex- 
periments, a diffuse magnetic background reminiscent of 
quantum disorder was observed while approaching Tn .— 
In NMR this magnetic background would have appeared 
as line broadening. Thus our experimentally observed 
narrow satellites with no further splitting do not indicate 
the presence of magnetic/structural disorder, which is in- 
consistent with the interpretations of /iSR (Ref. [29h and 
neutron scattering (Ref. |27| ) results measured on powder 
samples. At present the origin of this discrepancy is not 
clear to us. Further inelastic neutron scattering and /iSR 
experiments on single crystals are required for a better 
understanding of this issue. 

Our observations are very similar to ^'^As NMR mea- 
surements on (Ba,Sr,Ca)Fe2As2 in the stripe/columnar 
type magnetically ordered state^ii^i^ According to 
Kitagawa et al.r^ based on the crystal symme- 
try of (Ba,Sr,Ca)Fe2As2 with orthorhombic structure 
{Fmram), the internal field at the ^^As-site is explained 
by an off-diagonal transferred hyperfine field induced by 
four nearest-neighbor Fe moments whose direction is in 
the a5-plane. In the case of Neel type AFM order, the 
c-component of iJint at the As site is zero due to a per- 
fect cancellation of the off-diagonal hyperfine fields pro- 
duced by the four Fe moments when the spin moments 
are in the a6-plane. Only the stripe-type AFM order 
can produce a c-component of iJint at the As site. In the 



case of Pb2VO(P04)2, the crystal structure is monoclinic 
(P2i/a) with a lower symmetry. Thus the relationship 
between iJjnt and the magnetic structure is not trivial. In 
addition, the hyperfine field at the PI site contains both 
transferred and dipolar components. However, it may be 
possible to consider the present case based on the discus- 
sion for (Ba,Sr,Ca)Fe2As2 systems as a first approxima- 
tion because the local symmetry in the plane is close to 
four-fold axial symmetry corresponding to a square lat- 
tice formed by V ions with respect to the PI site. Since 
the magnetic easy axis is the &-axis and the direction of 
ifint is along the c-axis in the magnetic ordered state, it 
is concluded that magnetic order in Pb2VO(P04)2 has a 
columnar structure with spins in the a5-plane. Our in- 
terpretation in terms of a columnar magnetic structure 
is consistent with the neutron scattering experiments;^ 
Of course, as we pointed out, the four-fold symmetry of 
the square lattice is slightly broken due to the lower sym- 
metry of the crystal structure which makes the analysis 
more complicated. The small a- and 6-components of 
Hint might originate from the non-perfect cancellation of 
the hyperfine field due to the low symmetry crystal struc- 
ture. The observed splitting to four lines of the P2 site 
spectrum could also be related to such a complication. 

B. Dynamics 

The overall temperature dependent behaviors of 1/Ti 
for both P sites are alike. As shown in Fig. [H^b), 
1/ixTiT) above ~ 20 K is T- independent and increases 
slowly below 20 K where the system begins to show 
strong antiferromagnetic short-range correlations. In 
the same figure x(^) is also plotted to highlight the 
broad maximum corresponding to the short-range corre- 
lations. In general, \/{TiT) is expressed in terms of the 
dynamic susceptibility xa/(<?', "^o) per mole of electronic 
spins asiiiis^e 

A q 

where the sum is over wave vectors q within the first 
Brillouin zone, A{q) is the form factor of the hyper- 
fine interactions as a function of q in units of Oe/^B, 
and x^,^((f, tJo) is the imaginary part of the dynamic 
susceptibility at the nuclear Larmor frequency in 
units of iis/Oe. The uniform static molar suscepti- 
bihty X = XmC^jO) corresponds to the real compo- 
nent Xm('?*''^o) with q — and ujq — 0. Thus the 
temperature-independent behavior of l/(xTir) above 
20 K in Fig. [6Kb) demonstrates that the T-dependence 
of J2 I ^('?) P Xa/(9i'^o) scales to that of x- On the 

other hand, a slight increase of l/{xTiT) below 20 K in- 
dicates that J2 I ^(9) P Xm(9''^o) increases more than 
q 

X due to the growth of antiferromagnetic correlations 
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with decreasing T. There could be two possible sources 
for this growth of antiferromagnetic correlations. One is 
Neel type antiferromagnetic spin fluctuations character- 
ized by the wave vector q = (±7r/a, ±7r/6) and the other 
is columnar type AFM spin fluctuations with q — (±7r, 0) 
above Tn. The hyperfine form factor at the PI site can 
be written as |A((7)|^ = {2A[cos(ga;a/2) + cos{qyb / 2)]}'^ , 
if PI is located at the center of the perfect square lattice 
formed by V ions. In such a case, Neel type antiferro- 
magnetic spin fluctuations do not contribute to l/(xJiT) 

since is zero. However, in Pb2VO(P04)2, since 

PI is not located exactly at the center of a square lat- 
tice, A{q) should have a nonzero value. Therefore the 
enhancement of l/(xJir) could be due to the growth of 
Neel type AF spin fluctuations which are not completely 
filtered out at the PI site. Alternatively, since the AFM 
ordered state is columnar type as discussed before, if it 
has a finite correlation extended above Tn then [^(g)!^ 
would have a finite value which might also be responsi- 
ble for the enhancement of l/(xTiT). If large enough 
single crystals could be grown, the origin of the enhance- 
ment could be conclusively resolved from inelastic neu- 
tron scattering measurements of the magnetic fluctuation 
wavevector above Tn. 

A sharp peak in 1/Ti versus T is a direct indication 
of long-range magnetic ordering. However, the peak po- 
sition for the H = 0.65 T measurement is at 3.65 K 
while for 4.06 T it is enhanced to 3.93 K. We plot- 
ted these two points in the H vs. Tn phase diagram 
reported in Ref. [23 obtained from heat capacity mea- 
surements of Tn{H) [see inset of Fig. [Hb)], where it is 
clearly seen that Tn increases slightly with increasing H 
from to 7 T and then decreases with further increase of 
the field. This behavior is similar to that observed in FSL 
compound BaCdVO(P04)2 (Ref. [25l ) as well as in non- 
frustrated 2D antiferromagnet Mn(HCOO)2-2H20 [see 
p. 388 of Ref. Qualitatively, long-range magnetic 

ordering in low-dimensional and/or frustrated spin sys- 
tems is suppressed by quantum fluctuations. Magnetic 
field suppresses these quantum fluctuations, therefore Tn 
is slightly enhanced with increasing field. However, above 
10 T the field is strong enough to partially suppress the 
antiferromagnetic ordering, hence Tn is reduced. In the 
antiferromagnetic ordered state, 1/Ti is mainly driven 
by scattering of magnons off nuclear spins, leading to 
a power law temperature dependence4^'^ For T ^ A, 
where A is the gap in the spin-wave spectrum in tem- 
perature units, 1/Ti either follows a T^ behavior due 
to a two-magnon Raman process or a T^ behavior due 
to a three-magnon process, while for T ^ A, it follows 
an activated behavior 1/Ti cx T^exp(— A/T). As seen 
from Fig. [6la), our ^^P 1/Ti data in the lowest temper- 
ature region (1.5 K < T < 2.3 K) follow a T^ behavior 
rather than a T^ behavior suggesting that the relaxation 
is mainly governed by the two-magnon Raman process. 
The lack of activated behavior down to 1.5 K indicates 
that the upper limit of A is 1.5 K. Therefore, the up- 
per limit of the axial anisotropy [D w A^/(|Ji| -I- IJ2I)] 



(Ref. [H) would be about 0.15 K. 

The "^^P 1/Ti data and the hyperfine coupling con- 
stants determined above enable us to estimate the ex- 
change couplings in the system under investigation. Ac- 
cording to Moriya, the nuclear spin- lattice relaxation rate 
1/Ti at sufficiently high temperature in a system with 
exchange-coupled local moments is constant and can be 
expressed within the Gaussian approximation of the cor- 
relation function of the electronic spin as^^i^ 

X ^ ^ , (7) 



where uje = [max(|Ji|, \J2\)kB/h]y/2zS{S + is the 
Heisenberg exchange frequency, z is the number of 
nearest-neighbor spins of each V*^ ion, and z is the 
number of nearest-neighbor spins of the PI site. In 
order to get the coupling of the PI site to individual 
V^"*" spins, the total coupling is divided by z . Since 
our measurements were done for H \\ c, the hyperfine 
couplings in the a5-plane (i.e. A'^j: and A'f^j:) contribute 
to 1/Ti. Therefore, we have taken the rms average of 
couplings along the a- and 6-directions. The z in the 
numerator is due to the fact that the PI site feels the 
fluctuations arising from all nearest-neighbor V^^ spins. 
Using the relevant parameters {A^^j: ~ 1882 Oe//iB, 
Alf ~ 1251 Oe/fiB, z = 4, z' = 4, 5 = 1/2, and 
the high temperature (25 K to 100 K) relaxation rate 

of (^^^ ~ 127 s-i for the PI site) in Eq. 

the magnitude of the maximum exchange coupling con- 
stant was calculated to be max(|Ji|,|J2|) — 8 K. For 
a system with two different exchange mechanisms, the 
exchange frequency is mainly dominated by the larger 
exchange interaction. Therefore, in Pb2VO(P04)2, this 
value of J should represent the dominant exchange con- 
stant J2. Indeed, max(|Ji|,|J2|) — 8 K is close to the 
value J2 = (9.3 ± 0.6) K obtained from the K{T) analy- 
sis above. 



C. Critical effects 

The NMR spectrum measurements can give a pre- 
cise estimate of the temperature dependence of sublat- 
tice magnetization, which cannot be obtained from static 
magnetization measurements. The temperature depen- 
dence of Hint can potentially provide the critical expo- 
nent f3 reflecting the universality class/spin dimensional- 
ity of the spin system. In Pb2VO(P04)2, due to partial 
filtering of AF fluctuations above Tn, our 1/Ti data re- 
stricted us from doing a detailed analysis. The (3 values 
expected for different spin- and lattice-dimensionalities 
are listed in Table ^^L^L^M^ If we assume that the 
experimental P value of « 0.23 obtained above from 
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TABLE II: Critical exponent /3 expected for different univer- 
sality classes. The value of /3 listed for the 2D XY model is 
for a finite-size system. 







Ref. 


3D Heisenberg 


0.367, 0.33, 35 




3D XY 


0.345, 0.31, 0.33 


47^^ 


3D Ising 


0.326, 0.325, 0.31 


47^)53 


2D XY 


0.231 


54 


2D Ising 


1/8 


51,53 



the measurements on the single crystal is due to criti- 
cal fluctuations, then the (3 value is closer to the 2D XY 
than the 3D models although the ordering is suggested 
to be a 3D order from the powder neutron diffraction 
measurements.— We suggest that our data are not suf- 
ficiently close to Tn to reflect 3D critical behavior It 
is to be noted that the longitudinal /xSR relaxation rate 
above Tn in Pb2VO(P04)2 has been fitted well by this 
2D XY model. A tiny in-plane anisotropy can dominate 
the fluctuations if the dominant AF correlation length is 
largCf^ In the present case, x(^) above Tn is approxi- 
mately isotropic, but it is possible that there might be a 
small anisotropy which is not visible in the experimen- 
tal x(2^)- III fact, strong in-plane correlation may also 
enhance 3D coupling along the c-directioui^ Therefore 
to understand this issue, a precise estimation of in-plane 
spin anisotropy is required. 

It is useful to compare our data with data for other 
layered compounds. A similar reduced value of /3 close 
to the (3 for the 2D Ising model has been observed before 
in several layered compounds, e.g., K2NiF4, K2MnF4, 
K2C0F4, Rb2MnF4, and RbCoF4 (Ref. W). In these 
compounds the weak Ising-type anisotropy is about two 
orders of magnitude larger than the interplanar coupling. 
Therefore it is argued that the long-range ordering is pri- 
marily induced by the strong anisotropy crossover from 
2D Heisenberg to 2D Ising behavior, which precedes the 
2D Ising to 3D Ising crossover that should occur close to 
Tn. Since the onset of such 3D ordering is still governed 
by the 2D processes, the phase boundaries that are ob- 
served should reflect the underlying 2D character. When- 
ever there is a 3D ordering, one must have 3D correlations 
sufficiently close to Tn- For instance, in BaNi2(P04)2, an 
exponent (/3 = 0.33) close to 3D model was found near 
Tn and another exponent (/? = 0.23) close to that of the 
2D XY model was found a little bit away from Tn (see 
p. 281 of Ref. I47I ). However, in some cases this 3D criti- 
cal region is over such a narrow temperature range that 
it is not accessible experimentally. Such a scenario has 
been realized before in the frustrated antiferromagnetic 



square lattice compound Li2VOSi04. In this compound, 
neutron powder-diffraction experiments show that the or- 
dering is 3D and the antiferromagnetic a6-layers are cou- 
pled ferromagnetically along the c-axis.— In contrast, a 
reduced value of (3 close to the (3 value of the 2D XY 
model has been reported from NMR measurements and 
the authors suggested that the transition to the colum- 
nar phase might be driven by the XY anisotropy j^^i^^ 
Although our experimental observations below Tn are 
similar to those reported for Li2VOSi04, a proper ex- 
planation for this peculiar behavior is still lacking in all 
these vanadates. Therefore further experiments on high 
quality single crystalline materials are essential to ad- 
dress this issue. 



V. CONCLUSION 

We performed ''^P NMR measurements on 
Pb2VO(P04)2, which is a strongly frustrated 2D 
Heisenberg square lattice compound. The exchange 
couplings were estimated reliably from K{T) analysis 
to be Ji = (-5.4 ± 0.5) K and J2 = (9.3 ± 0.6) K. 
From the NMR spectral measurements, the ground 
state was detected to be columnar antiferromagnetic 
(CAF) type and is consistent with that expected from 
the phase diagram (Fig. [T]) for a = J2/J1 — —1.72. The 
hyperfine coupling and 1/Ti above Tn are consistent 
with the distortion of the V^"*" squares due to monoclinic 
crystal symmetry. We have measured the temperature 
dependence of the sublattice magnetization in the 
CAF ordered state from the line splitting. The critical 
exponent P estimated from the sublattice magnetization 
was close to the value predicted for a 2D XY model. We 
did not observe any signature of magnetic/structural 
disorder in contrast to the neutron scattering and /iSR 
experiments on powder samplesj ^'^'^^i^^ To gain more 
insight, further complementary experiments on single 
crystals are required. Due to the low energy scale 
of the exchange interactions, this compound is also 
a potential candidate for high-field experiments, and 
for high-pressure experiments where one can tune the 
exchange couplings. 
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